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Abstract 

We estimate that the LHC could set bounds BR{Z M*e=F) < 4.1 x 10"^ and BR{Z T^fJ."^) < 3.5 x lO"*^ 
(at 95% C.L.) with 20 fb~^ of data at 8 TeV. A similar sensitivity can be anticipated for Z — >■ r^e^, because we 

■ consider leptonic r decays such that Z — >■ t^^j,^ — >■ e*/i^+ invisible. These limits can be compared to the LEPl 

■ bounds of order 10^^ — )■ 10^". Such collider searches are sensitive to a flavour-changing effective Z coupling which 
^jqJ is energy dependent, so are complementary to bounds obtained from r £££ and /i — > eee. 
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I — ^. 1 Introduction 

^: 

Lepton flavour is conserved in the Standard Model (SM) of particle physics. However, the observation of neutrino 
^ I oscillations implies that there is Beyond-the-Standard-Model physics which changes lepton flavour. Lepton Flavour 
^ , Violation (LFV) among charged leptons is therefore a phenomenologically motivated place to search for elusive Beyond- 
^ ■ the-Standard-Model (BSM) physics. 

' The LHC is an electroweak-scale hadron collider, so can reasonably be expected to give information on electroweak 
fSJ ■ symmetry breaking, weakly interacting dark matter candidates (WIMPs) and new particles or interactions in the 
^ [ quark sector. In this paper, we are interested in the prospects of Lepton Flavour Violating processes at the LHC, 
\l . involving SM particles. LFV in BSM scenarios, such as supersymmetry or lieavy singlet neutrinos, has received more 
^ ; attention [1, 2, 3, 4]. 

■ The LHC produces many rs in the decays of W and Z bosons, and 6s (respectively ^ 10®, 3 x lO'', and 10^^ for 
• ' 10 fb^^ at 7 TeV [5]), so one could hope to be sensitive to LFV t decays such as t — > £j. LHCb obtained a bound 

^ ■ BR{t fififi) < 7.8 X IQ-^ with 1 Po^^ of data ^ [6]. However, for ATLAS and CMS, there do not seem to be enough 
^vq ^ high-pT rs, to allow an interesting sensitivity [5, 7]. We therefore explore LFV in the production of leptons — that is 
y—i • — in LFV decays of the Z ^ . The cross-section times branching ratio for pp Z ^ M^M^ is of order a nanobarn at the 
7 or 8 TeV LHC, so with 20 fb^^ of data at the end of 2012, the LHC wiU have 2 x 10^ BRiZ --5x10® Zs. 

This can be compared to the 1.7 x 10^ Zs at LEPl. Effective LFV Z couplings can also be constrained by the possible 
Z contribution to decays such as r — >■ /i/l/i and /i — !■ eee. In section 2, we review the complementarity of the collider 
and low energy searches. 

We study the decays Z — >■ ii^e"^ and Z — > r^/i^ (which results we extrapolate to Z r^e^). We focus on 
leptonic t decays, so the t^£^ final state is e^/i^ + two neutrinos. This is not without backgrounds, as discussed 
in section 3. The most dangerous is Z — )• r+r^ — > e*/^^ + 4 neutrinos, so section 2 also reviews some kinematic 
considerations, which suggest that this background could be reduced to an acceptable level. In section 4, we present 
our selection criteria, designed to maximise sensitivity to the signal. The bounds this analysis could obtain are in the 
final summary, which also discusses future prospects. 
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2 Theoretical review 

The three Z decays which violate lepton flavour conservation are Z e^/i^, Z — >■ /i^r^, and Z e^r^. We 
consider the detection of e^jX^ final states, arising either directly, Z — >■ e^/x^, or through leptonic tau decays: 
Z — > fi^T^ ^^e^i/D, and Z — > e^r^ — > e^fi^i'D. As the last two decay channels differ only by exchange 
of the electron and the muon, we only consider Z — > fi^r^ — > fi^e^i^D in this paper. We anticipate that the 
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^To be compared with the Belle bound in Table 2. 

■^LFV in the decay of Higgs produced at hadron colliders was discussed in [8], taking into account the constraints from low energy 
processes. An earlier study of the collider prospects was [9]. This subject has recently been revisited [10]. 
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sensitivity to Z ^ e^r^ should be equivalent, or better, because the SM backgrounds should be similar to those of 
Z — > fi^T^^ — > fi^e^uD, and the detection prospects for a low pt muon (from tau decay) are better than those of a 
low Pt electron. 



2.1 Flavour-changing Z couplings 

We suppose that new LFV particles arise at a scale M > mz, which allows to describe the flavour changing Z 
vertices with Effective Field Theory, via operators of dimension six and higher[ll]'^. The discussion applies to any 
flavour changing vertex, although we focus on the Zt~^ interaction to avoid flavour index sums. At dimension six, 
BuchmuUcr and Wyler[f 2] give the following SM gauge invariant operators which can contribute (with their hermitian 
conjugates) to Z ^ r^/U^: 

[H'<a'D„H][iyrir] , [H^D^H][I^j"er] , [H^D^H]-p^^T (1) 

I,a'Hafi,,TW'''^ , l^^HapcTB'^l^ (3) 

where SU(2) contractions are indicated in square brackets when there is more than one possbility, is the muon 
doublet, T the tau singlet, the hypercharge gauge boson, are the Pauli matrices, apa = §[7/3, 7a], B"^ = 
so on. We suppose that these operators are normalised such that, in the presence of electroweak 
symmetry breaking, the combination of operators on the first (or second, third) line of equation (1 -3) appear in the 
Lagrangian as the first (or second, third) line below: 

gzm%[Tl-io.{An,rPL + AR^rPR)Z°'T + h.c] (4) 

+2gz[ma{CL^rPL + CR^rPR)d^T + /l.C.]Z''« (5) 
+igzmr[lI(Jocp{Dl^rPL + -Dl^r^fl)T - /l.C.]Z«^ (6) 



where gz = VT+T^, ^"'^ = ^"^'^ - d^Z"", and the A,C,D coefiicients have mass dimension -2. We assume 
that these interactions are generated by a loop involving particles of mass M > mz, so A,C,D oc 5|/(167r^M^). 
Concretely, if Clij,t,Crht are taken real, and the Z is on-shell, then equation (5) and the SM Z couplings give 
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gv - gAl5 m%[CLnT + Cn^r - {Cli^t - Cn^rhb] 

m\[CL^T+CR^ir-{CL^tT-CR^r)lb\ 9v - 9a15 J V 

(7) 

where gv = —1/2 + 2s^r and gA = —1/2. 

The effective interactions relevant to LFV Z decays are those of equation (5). This is because the A and D 
coefficients are constrained by low energy data (see below), whereas the coefficients Cl,Cr are multiplied by the Z 
four-momentum-squared, which suppresses their contribution in low-energy tree^ processes. For instance, the effective 
flavour-changing coupling on the Z pole would be gzinziCLfirPL + Cr^tPr), whereas in r — >■ fijlfj, it would be 
^ gziTi^iCLfirPh + Cr^tPr)- Their contribution to Z decays is: 

BR{Z ^ r±M^) - 16m|(|C^^,|2 + \CR^r?) x BR{Z ^ m^M^) ~ 1-7 x 10"^^ (8) 

where — )• and the last estimate is for Cliit — Cr^t — ^^/(ISTr^M^). Recall that the effective operator formalism 
supposes m,z < M, so constraints ^ 10""'' are interesting. The current bounds on the Cs which can be extracted from 
the OPAL and DELPHf bounds on LFV Z decays [13, 14] are given in Table 1. 

Bounds on the C couplings can also be obtained from the contribution of a Z loop to decays such as r — )• /Lt7. 
To estimate a bound, we assume that M ^ mz, so the effective LFV vertices can be treated as contact interactions 
in the loop, and we suppose that at the scale M, New Physics induces C ~ g^/lQn^M^ and no A or D coefficients. 
Then, renormalisation group mixing from the scale M to mz can induce coefficients, which we estimate (for the 
case of T ^ /i7) to be 
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^Our discussion is patterned on [11]; additional references can be found in the references and citations thereof. 
^We thank Micheal Peskin for pointing out that this suppression may not occur in loops. 



process 


exp. limit 


16w'mlC/gl < 


BR{Z e±/i=F) (a) 
BR{Z^e^T^) (a) 
BR(Z II ' r^) (1)) 


1.7 X 10-« 

9.8 X K)-*^' 
1.2 X K)-""' 


0.32 
0.76 
0.84 



Table 1: Bounds on lepton flavour changing Z decays (a) from OPAL [13] and (b) from DELPHI [14]. The third 
colomn assumes the decay Z Z^Zj is mediated by CLa0 — Cr<x0 — fff /IStt^M^, and gives the resulting bound on 
m|/M^ (which is different for the various i^Zf )• 



The last colomn of Table 2 gives bounds on the -D^^s, which can translated into bounds on the C^^, 



67 decay imposes to|C^ ^ 10 °, which corresponds to BR{Z 



coefficients via 
< 10-10. The 



mz, so the loop calculation makes little sense, and the best bounds are the LEP limits of 



equation (9). The fi - 
r ^ £7 decays allow M 
table 1. 

We can make conflicting guesses for the "probable'' values of Cl and Cr. If the new particle mass scale M is 
beyond the reach of the LHC, then so are the LEV Z branching ratios (see equation 8). Also, as discussed below, the 
A and D coefficients are constrained by low energy rare decay data to be smaller than the LEPI bounds on the Cs. 
It is unclear to the authors how to build a model ^ with A, D <ti C. However, the Cs are physically distinct, so it 
is phenomenologically interesting to constrain them. Furthermore, the OPAL experiment observed a e+e" — >■ e+yu" 
event at LEP2 [15], which increases the interest of high centre-of-mass-energy searches for LFV. 

To obtain bounds from processes at scales ^ mz, the Z and its interactions (4-6) can be matched onto LFV four 
fermion operators. Current rare lepton decay bounds constrain the coefficients Al^tjAjiht more strictly ^ than the 
Z decay bounds of Table 1. The branching ratios mediated by the A coefficients are 
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BR{Z n^n'^) S 



\AL^.r? + \An^,r?) 
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BR{t — >■ nnn) 
BR{t^ 
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{2\AL^r? + \AR^r?) 



(10) 



so the rare lepton decay bounds of Table 2 put A^ and Ar at or beyond the sensitivity of the 7 or 8 TeV LHC. 

We neglect the D coefficients of equation (6), because their contribution to BR{Z t^/i^) is proportional to m^. 
This is different from the photon dipole operator 



where the factor is cancelled in BR{t fi'f) by a similar factor in the denominator. The SU{2) x U{1) gauge 

invariant operators of equation (3) induce both D"^ and , so D"^ ^ could be expected, in which case radiative 
lepton decays are more sensitive probes of dipole operators than Z decays. In Table 2, are given the current bounds 
normalised to g'^/16iT'^m%. 
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process 


exp. limit 


IQir'^m'zA/g'^ < 


16n'm%D/gl < 


BR{t -J> 666) 

BR{t — )■ /x/2/x) 

BR{li -> 666) 

BR{t 67) 
BR{t nj) 
BR{ii 67) 


2.7 X 10"** [17] 
2.1 X 10-** [17] 
1.0 X 10-12 [18] 

3.3 X 10-« [19] 

4.4 X 10"^ [19] 
2.4 X 10-12 [20] 


0.17 

0.15 
4.3 X 10-4 


1.3 X 10-3 

1.5 X 10-3 

4.6 X 10-6 



Table 2: Bounds from low energy Lepton Flavour Violating decays, on the A and D coefficients of equations (4) 
and (6). The third [fourth] colomn assumes the decay is mediated by ALfir — ^ii/^r — 5|/167r^M^ [-^l^t — 



g%/16Tr^M% and gives the bound on m^/M"^. If ef^ ~ gzD^, the bounds would apply to 
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could be suppressed with respect to C if the loop particles do not couple to the Higgs, but this does not suppress D couplings. 
®The operators inducing A coefficients also induce flavour-changing four fermion operators involving neutrinos. These are refered to as 
"Non-Standard neutrino Interactions", and can appear linearly in probabilities. Despite this potential enhanced sensitivity, the bounds 
[16] on IGir^ Am?2: / 9% generically ^ 1. 



2.2 Kinematics 



The LHC has more Zs than LEPl, but a less clean environment. So we here review some simple kinematics, for the 

decays Z t^£^, which suggest that the LHC backgrounds could be sufRcicntly reduced to probe couplings below 
the LEPl bounds. The analysis for Z — >• /x^e^ is relatively straightforward, since the /x^e^ should have the invariant 
mass of the Z, and is performed in section 4.2. 

In the case of Z ^ r^i"^, wc consider leptonic r decays, in particular — )• e^uP, so as to obtain opposite flavour 
charged leptons in the final state: 

Z -s- T^//^ {e^vV)ii^ . (11) 
The various backgrounds are studied in section 3; here we focus on the most problematic: 

Z t^t"^ {e^i^V){fifvV) , (12) 

which is refered to as "background" for the remainder of this section, and labelled Z + jets in Tables 3, 5 and 7. To 
estimate the number of such events, recall that BR{Z r^r'F) = 0.0337 and Bi?(r='= e^vv) = 0.178 [21]. 

With some simple approximations, one can show that the two-neutrino signal events of equation (11) may be 
kinematically distinguishable from the four-neutrino background events. First, since the rs are boosted by a factor of 
at least mz/^mr, it is not unreasonable to approximate as coUinear all the momenta of the r-daughters. For clarity, 
suppose that the t+ in the signal and background decays of equations (11) and (12) decays to e+. Then in both 
decays: 

Pt+ = Pe+ +Pu+Pu = Q!Pe+ • (13) 

Similarly, for the r~ of the background process, one can define p^- = (ip^- (where we neglect all final state masses). 
This gives: 

"^1 ^ = 2Q!Pe+ • p^- signal (14) 

m| — 2m^ = 2a/3pe+ " P^- background (15) 

As expected, the invariant mass 2p^+ -p^- of the charged final state leptons should be lower for the background events, 
and the pt of the signal muon will have a harder distribution than that of the r daughters (see Fig. 1). 

The second approximation is to neglect the transverse momentum pt,z of the Z, so that the pt of the Z-daughters 
balance. This implies that 

a\PT,e+ 1 = \Pt,ij.- I signal (16) 

o-\PT,e+\ = P\Pt,ix-\ background (17) 

To distinguish signal and background, we extract ai and a2 from our "data" using respectively equations (14) and (16), 
and use the difference 

Aa = a2 — cti (18) 

as a discriminating variable between background and signal (which peaks at zero). 

There is a small fraction of background events, with /3 — > 1, which constitute an irreducible background to our 
signal. These correspond to events where the /x carries (almost) all the momentum of its parent-r. They will occur 
when, in the reference frame of the decaying r, the /it carries close to the maximum allowed energy (m^ + mf^)/2mr, 
and when it is emitted in the direction of the boost of the t. If the charged lepton energies arc measured with a 
fractional uncertainty e, then we estimate that such irreducible events correspond to ~ of the background. For e ~ 
few %, this suggests that the background of equation (12) can be sufficiently reduced. 



3 Background and signal simulation 

The signal investigated in this paper consists of events with exactly one electron and one muon. At the LHC, the 
dominant SM process leading to e^/i^ dilepton final states is Z/j* — > t^t^ where both tau lepton decays Icptonically, 
one to a muon and neutrinos and the other to an electron and neutrinos. Two processes involving top quark(s) also 
lead to e^fjr^ final states: the top quark pair production channel pp ^ tt ^ M'^yM'~y and the single-top production 
channel pp — >■ Wt i^uM'^u with (.,£' = e, /U, r . Finally, gauge boson pair production can lead to e^/x^ states 

''The notation pp Wt includes both pp tW~ and pp — > tW~^ . 
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via: pp W+W- i+yl'-v and pp WZ/-f* I'vl+l' with = e, jx, r, pp WZ/Y qq'r+T', and 
pp Z/-i*Zh* ^ ///'/' with f^q,v, e, /i, T. 

The FEWZ [22] program was used to compute the inclusive pp — >■ Z/^* cross section at next-to- next-to-leading 
order (NNLO) at \/s = 8 TeV . The next-to-next-to-leading logarithm (NNLL) cross sections were taken from 
Ref. [23] for pp — >■ ti, and from Rcf. [24] for pp tW, with a top quark mass of 173.3 GcV. Finally, the diboson 
{pp — > WW, WZ, ZZ) cross sections were computed using mcfm [25] at next-to- leading order (NLO). These results 
are used to predict the total number of background events, which we then simulate at leading order. 

The simulation of those SM backgrounds and of the signal processes pp ^ Z ^ fi^r^ — > ji^e^vD and pp ^ Z ^ 
e^^jX^ was performed using pythia version 8.162 [26], except the simulation of the pp — >■ Wt — i^iyM'^v process which 
also used madgraphS version 1.4.6 [27]. In all Monte Carlo simulations, the cteq6l1 [28] parton density functions 
were used; initial and final state radiations and multi-parton interactions were also included. However, pile-up effects 
due to additional interactions were not taken into account in those simulations. 

All Monte Carlo events were passed through the delphes [29] program which provides a fast simulation of LHC 
detector response. The default configuration corresponding to the CMS detector [30] at the LHC was used. In the 
following, the pseudorapidity is defined as ?7 = — In [tan {0/2)] with the polar angle with respect to the beam direction; 
and (j) is the azimuthal angle. Isolated electrons and muons were identified by requiring no track with pt > 2 GeV in a 
cone of radius R = 0.5 around those lepton candidates, with AR = ^ Arf + Ac/)^. Standard Model backgrounds with 
non- isolated leptons were therefore not taken into account. Jets were reconstructed with the fast jet [31] program 
using the anti-kt algorithm with a distance parameter of 0.7. 

Table 3 contains the details of the Monte Carlo event production performed for the studies reported in this paper. 
As a total integrated luminosity of 20 fh~^ at ^/s = 8 TeV is assumed in the following analyses, the equivalent 
luminosities of these simulations are always higher than 200 fb~^ for each process to reduce statistical uncertainties. 





Processes 


ax BR (pb) 


Number of 








simulated events 


Z+]eis 


Z/7* ^ T±T=F l^vvl"^vv with 1,1' = e,ii {Mz/-y' > 20 GeV) 


237. 


47,500,000 


tt 


tt ->■ M+M'-u with I, e = e,ii,T 


23.6 


4,720,000 


single-top 


Wt P'vU'^v with e = e,fi,T 


2.35 


400,000 


diboson 


W+W- -^e+vi'-v with l,e = e,ii,T 


6.00 


1,200,000 




WZh* qq'T+T+ 


0.735 


140,000 




WZ/j* e^vl+l- with I, e = e,ii,T 


1.06 


220,000 




Zh*Zh* ^ ///'/' with / = g,i/,e, T 


8.26 


1,660,000 


signals 


Z -> T^iii^ e^vvii^ 




100,000 




Z n+e^ 




100,000 



Table 3: Standard model backgrounds in the search for Z decay violating lepton number at the LHC; the second 
column contains a x BR (in pb) for p-p collisions at ^/s = 8 TeV, and the third column shows the number of 
simulated events in this analysis which is in all cases greater than 10 times CaBR with £ = 20 fb~^. 



4 Sensitivity at the 8 TeV LHC 

For the signal investigated in this study, at least one lepton comes directly from the Z boson decay. As shown in 
Fig. 1, the distribution of this lepton pT has a clear peak around ^ 40 GeV, as expected for this jacobian distribution. 
At the LHC, a large fraction of SM backgrounds leading to e^/x^ final states comes W and Z bosons decaying to tau 
lepton(s) which decay leptonically (r — )■ iuv). In that case, the pr of the final lepton is strongly reduced due to the 
energy carried away by neutrinos (see Fig. 1). Requiring at least one high pT lepton, can significantly increase the 
signal to background ratio. This selection criterium also guarantees that events are accepted by single lepton trigger 
conditions of LHC experiments. In the following, 100% trigger efficiency is therefore assumed. 
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Figure 1: Lepton (electron or niuon) pr distribution at generator level at the 8 TeV LHC: for Icptons from pp — >■ 
Z/7* — )• (plain line) and for leptons from pp — > Z7* — > r+r^ — > i'^h'D Dv (dotted line), with £ = e or ji. 



4.1 The pp — >■ Z — > p,^T^ channel 

In the pp ^ Z ^ fi^T^ — > fi^e^Diy channel, the muon comes directly from the Z boson decay. The muon pT is 
required to be above 30 GeV and in the detector pseudorapidity acceptance I77I < 2.1. As the electron comes from the 
tau decay, the selection criterium on the electron pr has to be reduced to keep high signal efficiency. The electron 
PT is therefore required to be above 10 GeV and with |?7| < 2.5. To reduce the diboson background coming from 
events with more than 3 leptons, events are required to contain exactly 2 leptons. At this stage, a reducible diboson 
background comes from events where lepton(s) failed identification criteria leading to same-sign e^fj,^ final states. 
This background is therefore rejected by requiring two opposite-sign leptons. Double and single-top quark production 
are then strongly reduced by rejecting events which contain at least one jet with pr > 30 GeV and \r]\ < 2.5. To reduce 
the remaining diboson backgrounds, angular criteria in the transverse plane are used. For the signal, electrons and 
muons are almost back-to-back in the transverse plane, while the distribution of the A(j) angle between the muon and 
the electron is more flat for the diboson backgrounds. This A0(e,/i) angle is therefore required to be higher than 2.7 
radians. Another difference between backgrounds and signal is the A0 angle between the electron and the direction of 
the missing transverse energy (^t), which is peaked at for signal since electrons and missing transverse energy both 
come from a higly boosted tau lepton. This A(j>{e,]^T) angle is therefore required to be below 0.7 radians. Finally, the 
e/x invariant mass is required to be in the range [38 — 92] GeV, i.e. in the region around the broad peak expected for 
the signal. The numbers of events expected from SM backgrounds at each step of the event selection are reported in 
Table 4, together with the signal efficiency. 



Selection criteria 


backgrd. 


Signal efficiency (%) 


> 1 muon with pT > 30 GeV and jryj < 2.1 and 

> 1 electron with pT > 10 GeV and < 2.5 


48,181 ± 68 


9.4 ± 0.1 


exactly 2 leptons 

2 leptons with opposite charge 

jet veto: no jet with pT > 30 GeV and |?7| < 2.5 

A(j>{e,n) > 2.7 

Acj){e,^T) < 0.7 

38 GeV < Me,, < 92 GeV 


43,496 ± 65 
42,652 ± 64 
11,358 ± 34 
6,850 ± 26 
3,763 ± 19 
3,201 ± 18 


9.4 ± 0.1 
9.4 ± 0.1 

7.8 ± 0.1 

6.9 ± 0.1 
6.2 ± 0.1 
6.1 ± 0.1 



Table 4: Selection criteria for the Z — > r^fi^ search at the -y/s = 8 TeV LHC with C = 20 fh ^ with the total number 
of events expected from SM backgrounds and the signal efficiency (%); uncertainties are statistical only. 



The remaining SM backgrounds is 3201 events. As expected, 95% of this background comes from the Zj^* — >■ 
T^T^ -H- [i^er^vv background (see Tab. 5). The signal cfScicncy is 6.1 %. Normalized to the pp ^ Z cross section 
multiplied by BR{Z — > r^/x^) = 1.2 x 10~^, which corresponds to the current 95% C.L. limit from LEP, the number 
of events expected for the signal is 489. 



Processes 


Number of events 


Z+jets 


3,037 ± 17 


tt 


9 ± 1 


single-top 


7 ± 1 


diboson 


147 ± 4 


Total SM backgrounds 


3,201 ± 18 


Z — )• T*/i^ signal 


489 ± 6 



Table 5: After selection criteria of Table 4 for the Z t^ix^ search at the -v/s = 8 TeV LHC with £ = 20 /&~\ 
number of events expected for each SM background and for the signal normalized to cr{pp Z) x BR[Z r^/i^) 
with BR{Z — >• T^n^) = 1.2 x 10~^. These numbers correspond to Fig. 2. 



4.2 The pp^ Z ^ e^fj,^ channel 

In the pp ^ Z ^ e^/U^ channel, selection criteria are very close to those reported in previous section. The fact 
that this signal contains two high px leptons with a narrow invariant mass allows to reach a much better signal to 
background ratio. The muon and the electron pt are required to be above 30 GeV and their pseudorapidity |?7| < 2.1 
and < 2.5, respectively. Events are required to contain exactly two leptons, with opposite charge. The jet veto 
rejecting events with at least one jet with pr > 30 GeV and |?7| < 2.5 is then applied. Finally, the A0(e,/i) angle is 
required to be higher than 2.7 radians. The numbers of events expected from SM backgrounds at each step of the 
event selection are reported in Table 6, together with the signal efficiency. 



Selection criteria 


backgrd. 


Signal efficiency (%) 


> 1 muon with pT > 30 GeV and 77 < 2.1 and 

> 1 electron with pr > 30 GeV and Iry < 2.5 


36,489 ± 59 


18.5 ± 0.1 


exactly 2 leptons 

2 leptons with opposite charge 

jet veto: no jet with pt > 30 GeV and |?7| < 2.5 

A(l){e,n) > 2.7 


34,787 ± 58 
34,038 ± 58 
8,111 ± 28 
4,771 ± 22 


18.5 ± 0.1 
18.5 ± 0.1 
15.9 ± 0.1 
14.0 ± 0.1 



Table 6: Selection criteria for the Z e^/i^ search at the -^/s = 8 TeV LHC with £ = 20 /6 ^ with the total number 
of events expected from SM backgrounds and the signal efficiency (%); uncertainties are statistical only. 

The remaining SM backgrounds is 4771 events. At this stage the dominant background comes from Z/j* — > 
r^r^ — >■ ^.^e^uu events (see Tab. 7). And the signal efficiency is 14.0 %. Normalized to the pp ^ Z cross section 
multiplied by BR{Z e^/i^) = 1.7 x 10~®, which corresponds to the current 95% C.L. limit from LEP, the number 
of events expected for the signal is 158. In the mass region [87 — 95] around the Z peak, this signal would give 121 
events, whereas the total SM background represents 222 events, half of it coming r^r^ — )• n'^e^uu events 

and the other half from diboson events. 

4.3 Limits 

From the number of background events and the signal efficiencies reported in the previous section, expected 95% 
C.L. limits on the number of signal events can be computed and translated into limits on lepton flavour violating Z 

branching ratios. The modified frequentist CLs [32] method was used to compute those expected limits. The shapes of 
the Aa (defined in equation 18) and e^/U^ invariant mass distributions were used as final discriminant variables in the 
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Processes 


Number of events 


Z+jets 


2,756 ± 17 


tt 


199 ± 4 


single-top 


132 ± 4 


diboson 


1,683 ± 13 


Total SM backgrounds 


4,771 ± 22 


Z ^ e*/i^ signal 


158 ± 1 



Table 7: After selection criteria of the Z — )■ e^/x^ search at the -^/s = 8 TeV LHC with £ = 20 fb~^, number 
of events expected for each SM background and for the signal normalized to (j{pp Z) x BR{Z e^fJ-^) with 
BR{Z — >• e^/x^) = 1.7 X 10~®. These numbers correspond to Fig. 2. 



limit computation for the Z — > t ji^ and Z — ^ e channel, respectively. Those distributions are shown in Fig. 2. A 
3% systematic uncertainty was set on the SM background contributions and on the signal efficiencies. This uncertainty 
accounts for the experimental uncertainty on lepton identification [33], assuming that the uncertainty on luminosity 
measurement can be cancelled in our analysis using normalization to — >■ Z — ^ M"*"/^" ^^'^ PP ^ Z ^ e+e" events. 
This 3% systematic uncertainty was treated as fully correlated between signal and SM backgrounds. This procedure 
allows to set 95% C.L. expected limits of 3.5 x 10~^ and 4.1 x 10~^ on the Z boson branching ratios BR{Z — > r^/i^) 
and BR{Z — >■ e^/x^), respectively. Those expected limits are a factor ~4 better than the current experimental limits 
from LEP. 




A a e|x invariant mass (GeV) 

Figure 2: Search for flavour violating Z decay at the ^/s — 8 TeV LHC with C — 2Q fb^^: in the Z — > /i^r^ — > fi^e^vD, 
distribution of Aa after the final step of the event selection (left), and in the Z — >■ e^fi^ channel, distribution of the 
e^H^ invariant mass after the final step of the event selection (right). The dashed signal corresponding to the LEPl 
limits is superposed on the summed backgrounds. 



5 Summary 

The high number of Z bosons produced at the LHC allows to study rare Z boson decay predicted by BSM theories. 
Using 20 fb~^ of pp collisions at 8 TeV, we have shown that 95% C.L. limits on Z decay branching ratios to r^/i^ 
and e^fjr^ of 3.5 x 10~^ and 4.1 x 10^^, respectively, could be obtained. This represents a factor 4 improvements with 
respect to current best limits from LEP. The expected limit on BR{Z r^e^) should be equivalent to the one on 
BR{Z — >■ T^fj.^). As discussed in section 2, such collider searches are generically complementary to rare lepton decay 



searches, because they are sensitive to an effective vertex a s, which gives small contributions to tree level low energy 
processes. As shown in Eq. 8, reaching sensitivity to lower BRs would be interesting to test BSM models. However, 
an effective Z vertex (x s can also contribute in loops above mz', the limits on /x — >■ 67 are sufficiently restrictive to 
constrain this effect, implying that BR{Z M^e^) ^ 10~^°. 

In the analyses presented in this paper, basic and simple selection criteria were used to distinguish signals from 
SM backgrounds. More sophisticated analyses might further reduce backgrounds and increase sensitivities. In the 
limit computation, a 3% total systematic uncertainy on signal and background has been assigned, assuming some 
systematics, as the one on the integrated luminosity, could be cancelled in this measurement. 

It is interesting to speculate about the prospects of LFV Z decays at the post-2014 LHC. Our analysis is not 
immediately applicable, as the relative importance of the backgrounds will vary with the increased centre-of-mass 
energy. However, for the analysis described here, the expected limits reduce by factors of 2 or 3 with a luminosity 
of 200 fb~^, then further improvements are limited by the systematic uncertainties. Selection criteria which more 
effectively suppress the backgrounds, and reduced systematic uncertainties, would be required to profit from the full 
integrated luminosity of the future LHC. 
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